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The spontaneous emission of photons from optical cavities and from trapped atoms has been 
studied extensively in the framework of quantum optics. Theoretical predictions based on the rotat- 
ing wave approximation (RWA) are in general in very good agreement with experimental findings. 
However, current experiments aim at combining better and better cavities with large numbers of 
tightly confined atoms. Here we predict an energy concentrating mechanism in the behavior of such 
a composite quantum system which cannot be described by the RWA. Its result is the continuous 
leakage of photons through the cavity mirrors, even in the absence of external driving. We conclude 
with a discussion of the predicted phenomenon in the context of thermodynamics. 

PACS numbers: 42.50.Lc, 71.45.-d 



I. INTRODUCTION 

In quantum systems with purely unitary evolutions, 
the total energy is a conserved quantity, since the cor- 
responding Hamiltonian always commutes with itself. 
However, this argument does not apply to systems with 
Hamiltonians which are self-adjoint but not Hermitian 
[l[ and to open quantum systems [2]. The spontaneous 
emission of a photon is always related to the loss of energy 
from its source For example, in laser sideband cool- 

ing, a red-detuned laser field excites an electronic state of 
a strongly confined ion via the annihilation of a phonon 
from its motion. When followed by spontaneous photon 
emission, the phonon is permanently lost 6] . Contrarily 
to common believe, we show here that even un-excited 
and un-driven quantum systems might constantly leak 
energy into their environment. The origin of the pre- 
dicted effect are non-zero decay rates and the counter- 
rotating terms in the interaction Hamiltonian which are 
usually neglected as part of the RWA. 

Although the validity of this approximation has been 
questioned in the past 0, II]) it is commonly used to 
describe quantum optical systems. An exception is 
Hegerfeldt @ , who shows that the counter-rotating terms 
in the interaction between two atoms and the free radi- 
ation field can result in a small violation of Einstein's 
causality. Zheng et al. |10[ also avoided the RWA and 
predicted corrections to the spontaneous decay rate of a 
single atom at very short times. Recently, Werlang et 
al. [Hi] pointed out that it might be possible to obtain 
photons by simply placing an atom inside an optical cav- 
ity but no quantitative predictions have been made and 
no justification for the relevant master equation has been 
given. When avoiding the RWA, it has to be avoided ev- 
erywhere, also in the system-bath interaction. 

This paper contains a rigorous derivation of the quan- 
tum optical master equation for bosonic systems which 
uses only the Born and the dipole approximation. We ap- 



ply our results to individual quantum systems (trapped 
atoms and optical cavities) , and to a composite quantum 
system consisting of many atoms inside an optical res- 
onator. It is shown that for sufficiently large numbers 
of atoms inside the cavity, the stationary state photon 
emission rate can be as large as typical detector dark 
count rates. Its parameter dependence could be verified 
experimentally using currently available atom-cavity sys- 
tems Our calculations confirm the relevance of 
the effects predicted in [llj when the atom-cavity inter- 
actions are collectively enhanced [isl . [l6j . Similar energy 
concentrating effects might contribute significantly to the 
sudden heating in sonoluminescence experiments [17[ and 
are responsible for temperature limits in cooling experi- 
ments [ll [Hi ■ 

There are five sections in this paper. Section |lT] con- 
tains a rigorous derivation of the master equation of a 
single bosonic quantum system beyond the validity range 
of the RWA. Section Hn] calculates the corresponding sta- 
tionary photon emission rate. Assuming that the photon 
emission in the absence of external driving remains negli- 
gible in this case, we find that two of the constants in this 
master equation are approximately zero. Section HVl uses 
an analogously derived master equation to calculate the 
stationary state cavity photon emission rate for a com- 
posite quantum system consisting of many atoms inside 
an optical resonator. For feasible experimental parame- 
ters, we predict stationary state emission rates as high 
as 300 photons per second. A detailed discussion of our 
results can be found in Section fVl 



II. MASTER EQUATION OF A SINGLE 
QUANTUM SYSTEM BEYOND THE RWA 

Let us begin by studying individual quantum systems, 
like optical cavities and tightly confined atoms, with the 
ability to emit photons. Their Hamiltonian H in the 
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Schrodinger picture and in the Born and the dipole ap- 
proximation can be written a.s H — Hq + Hi^^ with [2 a ] 

Ho = hoj s+s" + '^hajk a,l^a,k\ , 

Hint = (^^^ akA+5kA flk^) s++h.c. (1) 

k,A 

and I^icaI = I^IcaI- The free radiation field consists of an 
infinite number of one-dimensional harmonic oscillators 
with wave vectors k, frequencies w^, polarisations A, an- 
nihilation operators OkA, and coupling constants gkA and 
^kA- In case of an optical cavity, w = Wc is the frequency 
of its field mode and the are the photon creation and 
annihilation boson operators and c: 



[c,ct] =1 , [c\c^] = 0= [c,c] , 



(2) 



In case of a large number of tightly confined two-level 
atoms with states |0) and |1), huj = Huq is the energy of 
the excited state |1) of a single atom and the are the 
collective raising and lowering operators 5*^ with 

[5+,^+] =0= , (3) 

as we shall see in the next paragraph. 

Suppose the atoms are confined in a region with linear 
dimensions that are much smaller than the wavelength 
of the emitted light. Then |k • (r^ — r^)! ^ 1 for most 
particle positions and Vj and for a wide range of wave 
vectors k. This implies that all particles experience ap- 
proximately the same g^x and ^kA- The Hamiltonian H 
remains therefore the same, if we replace s^s~ in Hq by 
CT3 and in Hint by cr^. Here cr^ and 173 are defined as 



N 
i=l 



0-3 



N 

E 



(4) 



with a3, = i (|1)..(1| - |0),,(0|), a+ = |1),,(0| and or ^ 
|0)ii(l| being the su(2) spin-like operators of atom i: 

, <jf] = ±af , [ai , a+] = -2a3. . (5) 

If the atoms are initially all in their ground state, they 
evolve under the action of the operators (|H) into the 
Dicke-symmetric states: 



lOp = [|0i0203...0Ar_aiV-i + llw-i+2...1) + 



+ |lll2...1iOi+iO,+2...0)]/ (^^^ 



1/2 



(6) 



which are the eigenstates of . The difference between 
excited and unexcited particles is counted by 173, since 
p(/|cr3||)p = / - In. For any I we have [H: 



lOp = ^N~{l-l)^^l \i- 



p I 



(7) 



This shows that and 0-3 are represented on |Z)p by the 
Holstein-Primakoff non-linear boson realization fl^ EB] 
CT+ = y/NS +As, g- ^ ^ AsS- with (73 = 5+5" - 
\N, As = ^1- S+S-/N, S+\l)p = VTTT\l + l)p, and 

'S' lOp = V^l^ ~ l)p foi' any I. The cr's still satisfy the 
su(2) algebra ([5]). However, for N I, ^ becomes 



(8) 



a± lOp = VNS^ \l) 



Consequently, in the large N limit, ([S|) contracts to the 
projective algebra e(2) [16] 



[53,5±]-±5±, [5-,^+]=l, 



(9) 



in terms of and S3 = as. This means, the oper- 
ators in the cavity case and in the many atom case are 
formally the same (cf. (II])-(l3l)). 

To derive the master equation of a single system, we 
assume that its state \if) at t = is known. Moreover, 
we notice that spontaneously emitted photons leave at 
a very high speed and cannot be reabsorbed. The free 
radiation field is hence initially in a state with only a 
negligible photon population in the optical regime 0- 
Denoting this state by \0), the (unnormalised) state 
vector of system and bath equals 



\0)\^",) = \0){0\ Ui{At,0) \OM 



(10) 



under the condition of no photon emission in (0, At). In 
the interaction picture with respect to Hq, this equation 
can be calculated using second oder perturbation theory, 
even when At ^ l/uj. Doing so, we find 

\ip1) = [l-As+s--B S-S+-C s+'^-D s-"^] \ip) (11) 

with 



At 



A = I At At' > gkA5kA e 

-^0 k.A 

I dt / di' Vg^;,5kAe-'(-+"^)(*-*') 



i(LJ~LJfc)(t-t') 



"'0 
At i-t 



k.A 



C = I At I At' VgkA5kAe'("-"'''*+'("+"'')*' 

'0 -^0 k,A 

At 

D = I At I Af 
'0 ^0 



Z.3kA5kA e-'(-+-'')*-'(— (12) 

k.A 



All four parameters could, in principle, be of first order 
in At due to the sum over the infinitely many modes of 
the free radiation field. 

In analogy to (jlOp . the (unnormalised) density matrix 
of the system in case of an emission equals 



Pi 



Tn. 



^a^akA Ui{At,0) pUl{At,0) 

k,A 



(13) 
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with p = \0){0\® p being the initial state of system and 
bath. Proceeding as above and using again second order 
perturbation theory, this yields 



Pi 



Bs^ps +C s ps +i)s+ps+. (14) 



The coefficients A, B, C, and D are obtained when tak- 
ing the complex conjugate of the coefficients A, B, C, 
and D in Eq. p2|) and extending the integration of the 
inner integral to At. To obtain relations between these 
coefficients, we decompose 



At j-At 

dt / dt'... = 

Ja 



At j-t 
dt / dt'. 

Jo 



At j-At 
dt / dt'. 

Jt 



(15) 

Substituting u = At — t and u' = At — t' in the second 
integral (which maps its area onto that of the first one) 
we find 



A = 2ReA , C = C* 
B = 2RcB , b = D* 



-2ii.jAt 



c. 



(16) 



Choosing the overall phase of C accordingly , the pa- 
rameters C, D, C, and D can hence be written as 



with 



C = D* = -f^c, C^D*^f*jc, (17) 



/eec''^'^* sm{ujAt)/uj (18) 



and with 7c being a real but not specified function of 
At. One can easily check that this notation is consistent 
with D = C* (cf. (US])) and with 1^. 

Averaging over the subensemble with and the 
subensemble without photon emission (cf. (fTTj) and (fT4|) ) 
at At hence yields the density matrix 



PI (At) = p 
1 
2 

-|-2Re-B s+ps" + 7c [/ s+ps+ + h.c] 



-^7C [(/ .s+^ 



-fB s s+) p + h.c] 
h.c.) p + h.c] +2Reyls"ps^ 



(19) 



In the following we return into the Schrodinger picture 
considering a master equation (28j 



^ = —n 



HcondP — pH^ 



nip) , 



' cond 

nip) = 7A s^ps^ +7B s+ps" -f 7c (s^ps^ +h.c.) , 

^/"cond = --?i[7A s+s" -I- 7b s~s+ -I- 7c (s'^^ -I- h.c.)] 

+riw s+s" (20) 

with 7A = 2ReA/At, 73 = 2ReB/At, and with uj be- 
ing the shifted bare transition frequency. Checking the 
validity of (1201) can be done easily by returning into the 
interaction picture and integrating pi(t) from zero to At. 
The result is indeed (III 



Before continuing, we remark that the master equa- 
tion (PU]) has been derived using second order perturba- 
tion theory. This means, it applies to quantum optical 
systems with the ability to spontaneously emit photons. 
The assumption of a Markovian bath has been avoided. 
Instead, we assumed rapidly repeated (absorbing) mea- 
surements whether or not a photon has been emitted 0- 
5j. These measurements constantly reset the free radi- 
ation field into \0) and make the system dynamics on 
the coarse grained time scale At with l/w ^ At ^ I/7 
automatically Markovian. Predictions based on this as- 
sumption have already been found in good agreement 
with actual experiments [2^ [2l| . 



III. PHOTON EMISSION FROM A SINGLE 
QUANTUM SYSTEM 

Let us continue by calculating the probability density 
= Tr(7?.(p)) for a photon emission of a system prepared 
in p. 



(7A s+s + 7b s ,s+ 



- 7C 



0+2 



')) . (21) 



Using ([2]) or ([3]), respectively, and considering the time 
evolution of the expectation values pi = (s+s"), ^1 = 
i(s~^ — s^'^'^), and ^2 = -I- s"*"^), we obtain a closed 
set of rate equations, 

Ml = -(7A - 7b) t*i + 7b , 

ii = -(7A - 7b) 6 + 2cj ^2 , 

6 = -(7A - 7b) 6 - 2^ a - 27c • (22) 

Setting these derivatives equal to zero, we find that 
the stationary photon emission rate of a single bosonic 
system (like an optical cavity or many tightly trapped 
atoms) is 



27A7B 27^ (7A - 7b) 



7A - 7b 4cj2 ^ _ ^g-)2 



(23) 



No photon emissions occur in the absence of external 
driving only when 



7B = 7c = , 



(24) 



as it is assumed almost everywhere in the literature 0- 
However, this assumption relies strongly on how 
the integrals in ([T^ are evaluated and whether relations 
like D{lu) = C{—Lu) are taken into account or not. 



IV. PHOTON EMISSION FROM A COMPOSITE 
QUANTUM SYSTEM 

Let us now have a look at a large number N of tightly 
confined atoms inside an optical cavity. The energy of 
this composite system is the sum of the free energy of 
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both subsystems, their interaction with the free radia- 
tion field, and the interaction between the atoms and the 
cavity field which changes ^ into 



Hq = h Uc c^c + hcuQ S^S + ^ hujk O-l^xCLkX , 

k,A 

-f^int = ^ HgkX OkA + .gkA OkA) + (9kA OkA 

k,A 

+9kA al^) S+ + y^hg^ (c + c^) 5-+ + h.c. (25) 

with 5c, ffkA, ffkA, gkA and gkA being coupling constants. 
For simplicity, the cavity photon states should be chosen 
such that gc becomes real. Proceeding as in the single 
system case, assuming the same properties of the bath, 
and returning into the Schrodinger picture we receive 
again the master equation ([20]) but with 



■Rip) 



K CpC 



:t + NT S+pS- , 



(26) 



Here cjc and ujq denote the bare atom and cavity frequen- 
cies, K is the cavity decay rate, and F is the decay rate 
of the excited state of a single atom. The crucial dif- 
ference to the usual Jaynes-Cummings model [22] is the 
presence of the cS" and the c'S^ term in ([25)) which 
vanish in the RWA. As we shall see below, these opera- 
tors result in a non-zero stationary state population in 
excited states and the continuous emission of photons, 
even without external driving. 

To calculate this rate, we take a conservative point of 
view and neglect 7b and 7c as suggested in (PH) . since 
this assures that no emissions occur in the absence of 
external driving in the single system case. Using ([2]), ([21), 
(I^Ul) . and (1251) we obtain again a closed set of the rate 
equations: 



m - 
m - 
m - 
m - 

6 = 
6 = 

with /i] 



'NgcTji - npLi , p2 = VNgcm - NTp2 , 
: 2\/]Vgc(l + 2^2 + + ^oVa + '^cVi - ^Cvi , 
: 2VNgc{l + 2pi + ^2) + i^oV4: + ^cV3 - ^Cm , 
-2VNgc{£,i + 6) - ^oVi - - ^Cva , 

2\/iVgc7?4 + 2wcC2 - , 

-2VNgcT]i - 2wcCi - , 

2VNg,ri4 + 2wo^4 - NT^s , 

-2VNg,ri2 - 2wo6 - ^r^4 (27) 



I, = {S+S-), m.2 = ± S+)ic T 



ct)),r,3.4=((5-T^+)(cTct)),^i=i(c2 
^3 ^ - 5+2), U ^ {S- 



42\ 
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5+2), and 
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FIG. 1: The cavity photon emission rate Ik as a function of 
time for A'^ = 10^ atoms inside the resonator obtained from 
a numerical solution of the respective rate equations. Here 
LOc =ujo = 384.2 ■ 10^^ s"^ {D2 line), = 6.1 • 10*s"\ F = 
1.9-10'^s'\ and k = 1.3-10^°s~\ as in Ref. [H]. After being 
initially empty, the cavity becomes populated — even in the 
absence of external driving. 



of these equations yields the cavity photon emission rate 



8C5c' + C'r + 4F i^o 



mCg^i^jQUc + 2CKr (22 (:d2) + 4kF (tig -^^y 

(28) 



( = K + NT. Combined with (I2T]) . the stationary state 



which applies for NT, \/Ngc, k <$i ujo,uJc- For example, 
the parameters of the recent cavity experiment with ^'^Rb 
[4] combined with N = 10* are expected to result in 
an as large as = 301 s~^ which can be detected 
experimentally (cf. Fig. [1]). 



V. CONCLUSIONS 

Let us now discuss the same topic in a more physical 
way. Among the reasons usually produced to justify the 
RWA, there is the one of the time scale, since the counter- 
rotating terms oscillate very rapidly so that their con- 
tribution to the dynamics remains in general negligible 
with respect to the resonating ones (23)]. Other authors 
apply the RWA in order to preserve quantum numbers 
and energy [23, ^255, since this approximation drops all 
the processes transferring energy between non-resonating 
modes. However, in Nature, there exist processes, where 
there is a redistribution of energy among different system 
degrees of freedom making possible some amounts of sys- 
tem self-organization. In particular, one could examine 
the possibility of concentrating the total energy of the 
system into a subset of degrees of freedom producing a 
decrease of its entropy, which in order to avoid a viola- 
tion of the second law of thermodynamics, would compel 
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the release of energy to the environment, thus keeping 
the free energy constant. This is possible, of course, only 
when the system is open. 

In this paper, we consequently examined the situa- 
tion occurring when the counter-rotating terms are not 
dropped. The avoidance of the RWA and the consequent 
more exact solution of the dynamics of the system under 
consideration gives rise to non-trivial consequences. As 
a matter of fact, the mechanism of concentration of en- 
ergy on a subset of degrees of freedom could help the un- 
derstanding of the hitherto rather mysterious processes 
of self-organization. The mathematical analysis we have 
done above shows indeed that in a quantum system a 
leakage of energy can occur among different degrees of 
freedom. This leakage is not necessarily triggered by an 
external pump of energy, but could be also triggered by 
the virtual photons coming from the quantum vacuum 
as, e.g., it occurs in the Casimir effect or in the Lamb 
shift. From the standpoint of the receiving system, the 
origin of the triggering energy is not important as far as 
the balance between the variations of energy and entropy 
is satisfied so to keep the free energy constant. In this re- 
spect, we recall that the ratio between these variations is 
just the temperature, as required by the thermodynamic 
definition: 

fcBT = ^ . (29) 

The interplay between the microscopic quantum dy- 
namics and the thermal properties certainly deserves fur- 
ther analysis, which, however, is out of the scope of the 
present paper. Here we are limiting ourselves to a spe- 
cific physical picture, which however, does not exclude 
other physical scenarios, as, for example, the conversion 
of energy from the thermal bath phonons to leaking pho- 



tons (as it occurs in the laser cooling mechanism). More 
appealing could be the occurrence of a dynamics, where 
a system is able to reach a state having a lower energy 
jumping over a separating barrier with the help of vir- 
tual photons coming from the vacuum. The possibility 
of realization of these scenarios needs of course further 
studies. First indications al ong the above lines can be 
found in the literature [l|, loMlllj. In all these examples, 
the system dynamics is irreversible. 

In summary, we derived the master equation for a sin- 
gle bosonic system (an optical cavity and a large number 
of tightly confined particles) without making any approx- 
imation other than the usual dipole and Born approxi- 
mation. We find that the effect of the counter-rotating 
terms in the interaction between a quantum optical sys- 
tem and its free radiation field might be annihilated 
by environment-induced measurements whether or not a 
photon has been emitted. Assuming that this is the case, 
we then show that these measurements cannot suppress 
the interaction between a large number of atoms and an 
optical cavity. The result is the continuous leakage of 
photons through the resonator mirrors, even in the ab- 
sence of external driving. For sufficiently many atoms, a 
relatively strong signal might be created. Its dependence 
on the system parameters can be verified experimentally 
using optical cavities like those described in [T2l-[l3|. We 
recognize that in order to better understand the physi- 
cal mechanisms responsible for the mathematical results 
presented in this paper, some more work is needed which 
is in our plans for future publications. 
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